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Abstract 

The Korteweg-de Vries equation lias a central place in a model for 
waves on shallow water and it is an example of the propagation of weakly 
dispersive and weakly nonlinear waves. Its history spans a period of about 
sixty years, starting with experiments of Scott Russell in 1834, followed by 
theoretical investigations of, among others, Lord Rayleigh and Boussinesq 
in 1871 and, finally, Korteweg and De Vries in 1895. 

In this essay we compare the work of Boussinesq and Korteweg-de 
Vries, stressing essential differences and some interesting connections. Al- 
though there exist a number of articles, reviewing the origin and birth of 
the Korteweg-de Vries equations, connections and differences, not gener- 
ally known, are reported. 

A.M.S. Classification: Primary 01-02, 01A55; Secondary 76-03, 76B25, 
35Q53. 
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1 Introduction 



It was in the "interest of Higher Truth" that professor Martin Kruskal, at 
the conference in commemoration of the centennial of the pubhcation of the 
Korteweg-de Vries paper in the Philosophical Magazine claimed that "he, 
together with professor Norman Zahusky, was the person, who more than any- 
one else, resuscitated the Korteweg-de Vries equation after its long period of, if 
not oblivion, at least neglect", Indeed it is well-known that in the follow-up 
of their 1965 paper in the Physical Review Letters "Interactions of Solitons 
in a coUisionless plasma and the recurrence of initial states" , a real explosion 
of research on this and related equations appeared in the journals. Many de- 
velopments in several fields of pure and applied mathematics, physics, chem- 
istry, biology and engineering followed. Restricting to mathematics we mention 
analysis, integrability of nonlinear systems, Lie-algebra's, differential geometry, 
quantum and statistical mechanics, PIEIE]. 

In this essay we direct our attention to the origin of the Korteweg-de Vries 
equation and its birth which has been a long process and spanned a period of 
about sixty years, beginning with the experiments of Scott-Russell in 1834 
the investigations of Boussinesq and Rayleigh around 1870 [7-11] and finally 
ending with the article by Korteweg and De Vries in 1895 

In simplified form the Korteweg-de Vries equation reads 

du „ du d'^u 

and it is the result of research concerning long waves in shallow water; x and t 
denote position and time and u = u{x,t) the wave surface. 

Nowadays, it is hard to understand for mathematicians not specialized in 
fluid mechanics, that a subject as this could raise such a wide spread interest. 
However, this was not the case in the nineteenth century when the study of 
water waves was of vital interest for applications in naval architecture and for 
the knowledge of tides and floods. Notably in England and France much research 
was spent on the study of water waves of several kinds, in England by, among 
others, Scott-Russell, Airy, Stokes, McCowan, Lord Rayleigh and Lamb and in 
France by Lagrange, Clapeyron, Bazin, St. Venant and Boussinesq. 

In some treatises and textbooks on "soliton theory" a short survey of the 
early history is presented but apparently it has not been the intention of 

the authors to dwell extensively upon the considerations and the mathematical 
analysis of those present at the cradle of the equation that became later known 
as the Korteweg-de Vries equation. Nevertheless, there are some review papers 
where more specific attention has been given to investigations related to the 
Korteweg-de Vries equation. We mention in particular the reviews by BuUough 
[Tl|. Bullough and Caudry Miles ^Hl and the recent impressive extensive 
article by Darrigol [17] and a letter in the Notices of the A. M.S. by Pego |18| . 

In these articles the work by Boussinesq on the one side and that of Korteweg 
and Dc Vries on the other side have been discussed. Studying these papers, 
the present author became aware of some inaccuracies regarding the relevance 
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and the significance of the work by Korteweg and De Vries, maybe even a 
slight animosity over the priority of the discovery of the equation. For example, 
Miles, Darrigol and Pego suggest that Korteweg and De Vries were presumably 
unaware of the work by Boussinesq. This is to be doubted because in their 
article reference has been made to the Comptes Rendus papers by Boussinesq 
[7] and St. Venant Besides this, the historian B. Willink [20] has presented 
the author with a handwritten copy of De Vries, containing an excerpt of the 
paper by St. Venant ^\ and there appears a clear reference to the "Essai sur 
la theorie des eaux courantes" ^HIj which proves that De Vries was certainly 
aware of the existence of Boussinesq's research. Pego writes "It is not clear 
why Korteweg and De Vries thought the permanence of the solitary wave still 
controversial in 1895" This is in contrast with the introduction of the 

KdV article, where it is stated "They (Lord Rayleigh and McCowan) are as it 
seems to us, inclined to the opinion that the wave is only stationary to a certain 
approximation. It is the desire to settle this question definitively which has led 
us into somewhat tedious calculations, which are to be found at the end of our 
paper" , HI. 

It is evident that Korteweg and De Vries, wanting to check the theory of 
long waves in shallow water, use an independent approach. It is our intention 
to illustrate this in the next sections, pointing out not only differences but also 
close connections in both theories. Here we already give some examples. Boussi- 
nesq used a fixed coordinate system and Korteweg and De Vries a coordinate 
system moving with the wave. The central equations in Boussinesq's analysis 
are the continuity equation and an expression for the wave velocity [H], whereas 
the Korteweg-de Vries equation is the central equation to which Korteweg and 
De Vries frequently revert in the course of their paper ^ . A simple substitu- 
tion of the wave velocity into the continuity equation yields immediately the 
Korteweg-de Vries equation in its full glory. However, Boussinesq did not do 
this, otherwise it may well be that the history of the long stationary wave had 
taken a different course. Pego ^Sl pointed out that the Korteweg-de Vries equa- 
tion appeared already in a footnote on page 360 of Boussinesq's 680 pages vast 
volume "Memoir sur la theorie des eaux courantes" ^U], that appeared in 1877, 
well before the publication of the Korteweg-de Vries paper in 1895. However, 
this footnote on the Korteweg-de Vries equation and also Boussinesq's expres- 
sion for the wave velocity are only valid when the wave vanishes at infinity, 
while this is not necessary in the theory of Korteweg and De Vries. Therefore, 
Boussinesq uses another approach for treating steady periodic waves than Kor- 
teweg and De Vries, who presented a unified treatment for steady waves, not 
only vanishing at infinity but for waves being periodic as well. It is not only the 
equation , but also its applicability that is important. It seems that this is not 
always sufficiently realized or even mentioned in the literature. Darrigol spends 
in his essay only one page to the Korteweg-de Vries equation under the heading 
"The so-called Korteweg-de Vries equation" |17| . It is only a whim of Tyche, 
the daughter of Zeus and the personification of fate, that Zabusky and Kruskal 
attributed the names of Korteweg and Dc Vries to our equation and not that of 
Boussinesq, who merits of course the token of priority. 



3 



In the following account we present a review of the work by Boussinesq and 
Korteweg and De Vries; as to Boussinesq, most of our attention is directed to his 
long article in the Journal des Mathematiques Pures et Appliquees , which is 
more accessible than his vast memoir ^U). The author, not a historian, is well 
aware that he may have overlooked or deleted important facts, but nevertheless 
he hopes that this study may disclose some generally unknown aspects of the 
early history of the Korteweg-de Vries equation. 

2 Scott Russell's experiments 

The story of the discovery of the "Wave of Translation" by John Scott Russell 
in 1834, has been recorded in many books concerning "Soliton" theory, the 
more so because Scott Russell's account is fascinating and even full of emotion, 
hardly expected in a scientific paper. Therefore, our account will be rather short 
and the interested reader is referred to his "Report on Waves" [SI HI El ■ 
was in the year of 1834 that the Scottish naval architect followed on horseback a 
towboat, pulled by a pair of horses along the Union Canal, connecting Edinburgh 
and Glasgow. However, the boat was suddenly stopped in its speed - presumably 
by some obstacle - but not the mass of water, which it had put in motion. Our 
engineer perceived a very peculiar phenomenon: a nice round and smooth wave 
- a well defined heap of water - loosened itself from the stern and moved off 
in forward direction without changing its form with a speed of about eight 
miles an hour and about thirty feet long and one or two feet in height. He 
followed the wave on his horse and after a chase of one or two miles he lost 
the heap of water in the windings of the channel |SJ . Many a physicist would 
not be inclined to analyze this phenomenon and leave it as it is, not so Scott 
Russell discovering something very peculiar in a seemingly ordinary event. He 
designed experiments generating long waves in long shallow basins filled with 
a layer of water and he investigated the phenomenon he had observed. He 
studied the form of the waves, their speed of propagation and stability, clearly 
perceptible in progressing positive waves, but not in progressing negative waves. 
A schematic view of these experiments is shown in figure 1, which is adapted 
from Remoissenet 




Figure 1: Scott Russell's experiments 
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For an extensive historical study the reader may consult the papers by Bul- 
lough and Darrigol ^Tj. As mentioned in the introduction, there existed in 
England and France a rich tradition in the mathematical description of hydro- 
dynamic phenomena such as wave motions in fluids. Scott Russell challenged 
the mathematical community to prove theoretically the existence of his soli- 
tary wave and "to give an a priori demonstration a posteriori" i.e. to show the 
possible existence of a stable solitary wave propagating without change of form. 

It is not unusual that new discoveries or new ideas encounter resistance 
from established convictions. We take from Rayleigh's paper "On Waves" in 
the Philosophical Magazine [11, pp 257-279, 1876] the following quotations. 
Airy, an authority on the subject, writes in his treatise on "Tides and Waves" 
|22|: "We are not disposed to recognize this wave (discovered by Scott Russell) 
as deserving the epithets "great" or "primary", and we conceive that ever since 
it was known that the theory of shallow waves of great length was contained in 
the equation — g^^^ the theory of the solitary wave has been perfectly 
well known" . Further "Some experiments were made by Mr. Russell on what 
he calls a negative wave. But (we know not why) he appears not to have been 
satisfied with these experiments and had omitted them in his abstract. All of 
the theorems of our IVth section, without exception, apply to these as well as 
to positive waves, the sign of the coefficient only being changed". Probably it 
was also Airy who expressed for the first time as his opinion that long waves 
in a canal with rectangular cross section must necessarily change their form as 
they advance, becoming steeper in front and less steep behind and in this he 
was supported by Lamb and Busset P |221 ■ Stokes believed that the only 
permanent wave should be basically sinusoidal, but later on he admitted that 
he had made a mistake, see also our section 7. 

On the other hand he writes |24| : "R is the opinion of Mr. Russell that the 
solitary wave is a phenomenon "sui generis", in no wise deriving its character 
from the circumstance of the generation of the wave. His experiments seem to 
render this conclusion probable. Should it be correct, the analytical character of 
the solitary wave remains to be discussed". 

The "a priori demonstration a posteriori" asked for by Scott Russell was 
finally given, first by J. Boussinesq in 1871 [7-10], some time later in 1876, by 
Lord Rayleigh [llj and in order to remove still existing doubts over the existence 
of the solitary wave by G. de Vries [23 and by D.J. Korteweg and G. de Vries 
in 1895 n. 

In the next section we present first a concise account of the contribution by 
Rayleigh, since it is short and leads directly to the heart of the matter. Moreover, 
this paper has been of great influence on the research of Korteweg and De Vries. 
Consecutively, we discuss in the other sections the investigations of Boussinesq 
and Korteweg-de Vries and we finish with some concluding remarks. 
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3 Rayleigh's Solution 

Be given an incompressible irrotational flow in a canal with a constant rectan- 
gular cross-section, fig. 2. 




Figure 2: Wavesurface 

The coordinates of a fluid particle are given by the coordinates x and y, 
the undisturbed depth of the canal by the constant H and the wave surface by 
H + h{x,t). Another essential assumption is that the wave length is large in 
comparison with H. 

As has been mentioned in the preceding section, Airy had already remarked 
that the theory of shallow waves of great length is contained in the equation 

(where k = H and g the constant of gravity) . The wave velocity is ^/gK , a result 
generally known since Lagrange in 1786. However, this result is only valid as a 
first order approximation, where h/H may be neglected. 

Rayleigh remarked that for this value of the wave velocity the so-called 
free surface condition (equilibrium of pressure) is only satisfied whenever the 
ratio h/H may be neglected, but if this is not the case it is impossible to 
have a wave in still water with velocity -^/gw and at the same time propagating 
without change of form. In order to cure this discrepancy with Scott Russell's 
experimental results, he proposes to look for a more accurate approximation 
of the wave velocity (dl, pp 252-253). Rayleigh assumes the existence of a 
stationary wave, vanishing at infinity, and by adding to the fluid a yet unknown 
constant basic velocity equal and opposite to that of the wave, he may omit the 
dependence on time. Since the flow is free of rotation, and incompressible, there 
exist a velocity potential (f) and a stream function ip, both satisfying Laplace's 
equation. The horizontal and vertical velocity components are given by 

d(j) dip d<j) dip 

u = — = "H- and — = - — 

ax ay ay ax 

and a series expansion gives 

u = 1^ = + (3-1) 

,3 



i^ = -ym + |/<^' (-)-... (3.2) 
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This expansion is justified because f{x) is, due to the large wave lengths, a 
slowly varying function of x. 

Integration yields the stream function 

V' = 2//(^)-|^/"(x) + |^/(^)(x) (3.3) 

constant along stream lines and hence also along the wave surface y = H + h(x). 
Let p be the pressure just below the wave surface, then we have the relation 

-2^—— = 2g{H + h) + ul + v^^:=Q (3.4) 

where p is the density, g the constant of gravity and C an integration constant; 
the suffix s denotes that the values of u and v are taken at the surface of the 
wave. 

To satisfy in higher approximation the free surface condition, — p constant 
— , Rayleigh investigates how far it is possible to make uj constant by varying 
h{x) as function of x. Using + = u'^{l + (^)^} and eliminating the 
unknown function / with the aid of l|3.1() and (|3.3|l he obtains after a tedious 
calculation a differential equation for the wave form y = H + h{x), viz. 

where uq is the still unknown constant basic velocity that has been added to 
the flow and C is again an integration constant, (pJJ, pp 258-259). The cubic 
expression at the right hand side vanishes for y = H with a; = oo and for 
y = H + ho with hg the crest of the wave. 
Elimination of C yields 

uo = VgiH + ho) « ^+^ho^ (3.6) 
which is also the wave velocity, and the equation (|3.5|l reduces to 

with the solution 



hix)^hosech'\J^xj. (3.7) 

This formula represents the "heap of water" (the "Great Wave") with the 
right wave velocity 1)3. 6|l , as experimentally observed by Scott Russell and so he 
was finally vindicated after fourty years of much discussion in England. 

Rayleigh finishes his article with the remark: "/ have lately seen a memoir 
by M. Boussinesq, Comptes Rendus, Vol. LXXII, in which is contained a theory 
of the solitary wave very similar to that of this paper. So as far as our results 
are common, the credit of priority belongs of course to M. Boussinesq". 
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4 The Equations of Boussinesq 



From Miles we learn in his essay ^B] that Boussinesq (1842-1929) received his 
doctorate from the Faculte des Sciences, Paris, in 1867, occupied chairs at Lille 
from 1873 to 1885, and at the Sorbonne from 1885 to 1896. He made significant 
contributions to hydrodynamics and the theories of elasticity, light and heat. He 
wrote several papers on nonlinear dispersive waves |7| El El ^tnd a voluminous 
"memoire" , entitled "Essai sur la theorie des eaux courantes" , presented to the 
Academic des Sciences in 1877, Vol. XXIII (JjQ;, pp 1-680). 

The most accessible publication is his article in the Journal de Mathema- 
tiques Pures et Appliquees in 1872 0. It has the verbose title "Theorie des 
ondes et des remous qui se propagent le long d'un canal rectangulaire horizontal, 
en communiquant au liquide continu dans ce canal des vitesses sensiblement 
pareilles de la surface au fond". This paper subsumes the short Comptes Rendus 
[3 IH] , whereas the monograph ^U] gives also much information less relevant for 
this exposition concerning the "Great Wave" . 

He considers, in the same way as Rayleigh, long waves in a shallow canal 
with rectangular cross section; the fluid is supposed incompressible and rotation 
free, while friction, also along the boundaries of the canal, is neglected. Distinct 
from Rayleigh's article, Boussinesq introduces also a time variable, essential for 
the description of a dynamic phenomenon, and the coordinates of a fluid particle 
at time t are denoted by {x,y) — {x{t),y{t)), see fig. 2. 

Be p the pressure in the fluid, p its density and (w, v) the velocity vector. 
The height of the water in equilibrium is again denoted by the constant H and 
the wave surface by the function y = H + h{x, t). The wave length is supposed 
to be large and the amplitude h of the wave small in comparison with H and 
vanishing for a; — > ±oo. 

Boussinesq's exposition starts along the same lines as in the theory of Ray- 
leigh. The main ingredients are a series development into powers of y, similar 
as in (EU - 

= + (4.1) 

with the as yet unknown function / = f{x, t) and valid for Q < y < H -\- h{x, t). 
Further, he uses the free surface condition 

gh+\{ul+vl) + ^+x{t)^0 (4.2) 

where x(0 is an arbitrary function and where the suffix s refers to the wave 
surface. Under the assumption that the potential 4> and its derivatives with 
respect to x, y and t vanish for x — > ±oo, the function x(i) may be omitted. 
We shall see that this assumption is very essential in Boussinesq's theory and 
it is used again and again in his paper, see also section 7. 

A second boundary condition follows from the kinematic equation 

dh dh dh 
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Substitution of the series expansions of the potential and the velocity com- 
ponents into (|4.2() and 1)4. 3|l results into two equations, containing h{x,t) and 
f{x,t). Elimination of f{x,t) gives in a first approximation, with h{x,t) small 
in comparison with H , the wave equation of Lagrange: 



A second higher approximation yields the well known equation of Boussinesq 

2H ^ Tdx^ 



dx^ 



(4.4) 



This equation may be simplified by restricting the theory to waves propagating 
into only one direction, say the positive x-axis. Denoting the wave velocity by 
w(x,t) and using the conservation of mass 



dh ( ii) /j^^ 
dt dx dt dx dx 



(4.5) 
we obtain after substitution into (|4.4() and integration with respect to x 



d_ 

dt 



(ujh) 



2ll 



d^h 
~d^' 







(4.6) 



To make progress it is desirable to have an explicit expression for uj{x,t), 
because substitution into H4.5() or (|4.6() gives a differential equation for the wave 
surface h(x,t). To this end Boussinesq introduces without a clear motivation 
the function 



2 



d^h 
3 dx'^ 



oT7 + ^7^ ■ (4-7) 



Differentiation of this expression with respect to t and replacing ^ by 
—\JgH-^, which does not disturb the order of approximation in the second 
term, yields 



Substitution of H4.6(l and (|4.5|l gives = a 



3 dx^ 
from which it follows 



that = for a wave propagating in the positive x-direction. Hence, Boussinesq 
obtains from H4.7(l the important result 



u;{x,t) = y/gH + y/ gH 



I Iff ^ 'Ghd^J 



(4.8) 



From 1)4. 5|) and 1)4. 8|) one may obtain the differential equations 



dh 
I— 
dt 



dh 
'dt 



- U! 



dh 

dx 



, duj 
dx 



d_ 

dx 



(h^Lo) + 2ujh 



dh 
dx 
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1 



and 

dh _ 1 Ay 1 (9 
H ~ ~4:\jHhdx 
or after passing to the new variable h dx 

dh 1 Ag" d 
'dt ~ IVHd^ 



3 h \dx h dx 
—da 



2 
3 



(4.9) 



(4.10) 



It appears from H4.8|) that the wave velocity differs from point to point at 
the wave surface and so it is expected that the wave should change its form 
during its course, which is one of the main issues in pursuing Scott Russell's 
experiments. However, the wave will only be stationary whenever a;(x,t) is 
constant. It is appropriate to make here some comments: 

1. The introduction of the function ip in H4.7|l is not a priori clear and 
motivated. With the aid of (|4.5|) wc may write instead of (|4.t)|) 



0. 



Replacing ^ by —y/gH and integrating with respect to x, we get the result 
H4.8II : remember h(x,t) and its derivatives vanish for x ±oo. 

2. Another proof of 14.81) was given many years later in 1885 by St. Venant 
|19|. It may be that he was not satisfied with Boussinesq's derivation. He 
applied another approach, using the mean value of the horizontal component of 
the velocity vector 



U{x,t) = 



1 



H 



H+h{x,t) 



u dy 



h{x,t) Jo 

3. The equations (|4.9|l and (|4.1UI) contain implicitly the wave velocity cu. 
This dependence on uj can simply be eliminated by the substitution of (|4.8|l 
into (|4.5|) . If Boussinesq had carried out this small operation he had obtained 
the Korteweg-de Vries equation long "avant la lettre" , viz. 



dh 
'dt 



^'-K'-Hh+h^ 
H2dx \3 2 



d^h 
~9"d^ 



= 



(4.11) 



valid for waves vanishing at infinity. This equation does not differ essentially 
from the Korteweg-de Vries equation as presented in the Korteweg-de Vries 
paper in the Philosophical Magazine; the coordinates in H4.11|l refer to a fixed 
(a;, t) frame, whereas Korteweg and De Vries used a moving frame, see also next 
section. 

4. As mentioned in the Introduction of this paper, R. Pego jJHI and also 
O. Darrigol (JT], pp 47) have discovered in a footnote on page 360 of the "Essai 
sur la theorie des eaux courantes" , that Boussinesq had found already in 1876 
the Korteweg-de Vries equation, mind without recourse to the expression (|4.8|l 
for the wave velocity ^Ol- In fact he used instead of H4.7(l the function 

> . ^ dh r-Frdh 1 ^9 f3h^ H"^ d'^h^ 



dx 



dx V 2H 



3 dx^J 



chp_ 
dx 
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and using ip = he gets equation H4.11|l . 

Consecutively, the wave velocity Lu{x,t) may be determined by the integra- 
tion of (|4.5|l . i.e. 

5. The bidirectional Boussinesq equation l|4.4|l can be factorized as 

\jt - n^Tx (a^'^ + 2^^ + j h « ^^'^^^ 

from which it immediately follows that the unidirectional Korteweg-de Vries 
equation 1)4.11(1 is contained in Boussinesq's equation 1)4. 4|l . 



5 The Korteweg-de Vries Equation 

We start with a few biographical data of Korteweg and De Vries [201 ■ Diederik 
Johannes Korteweg (1848-1941) received in 1878 his doctorate at the University 
of Amsterdam, after defending his thesis on the propagation of waves in elastic 
tubes. His supervisor was J.D. van dcr Waals, renowned for his equation of 
state and the continuity of the gas and fluid phases. Korteweg occupied the 
chair of mathematics, mechanics and astronomy at the University of Amster- 
dam from 1881 to 1918; he published several papers on mathematics, classical 
mechanics, fluid mechanics and thermodynamics. We mention in particular his 
investigations on the properties of "folded" surfaces in the neighbourhood of 
singular points, work related to that of Van der Waals Another scientiflc 
achievement is the edition of the "Oevres Completes" of Christiaan Huygens 
and Korteweg was the principal leader of this project in the period 1911-1927. 

He inspired many young mathematicians who wrote their thesis under his 
supervision, among others Gustav de Vries and the famous L.E.J. Brouwer. 
Korteweg had a great influence on academic life in the Netherlands as appears 
from his leadership in several academic institutions. The thesis of Gustav de 
Vries, entitled "Bijdrage tot de Kennis der Lange Golven" was published 
in 1894 and the paper in the Philosophical Magazine of 1895 is an excerption 
of this thesis. De Vries published papers on cyclones in 1896 and 1897 and two 
papers "Calculus Rationum" in the proceedings of the Royal Dutch Academy of 
Arts and Sciences in 1912. He teached mathematics at a gymnasium in Alkmaar 
and at a secondary school in Haarlem. 

The author received from the grandsons of De Vries a copy of the scien- 
tific correspondence between Korteweg and De Vries. From this we know that 
Korteweg advised De Vries to study Rayleigh's method of the series expansion 
which has been explained in section 3 of this paper. He also suggested to include 
capillarity and to investigate long periodic waves. 

Korteweg and De Vries start their article with the time dependent Rayleigh 
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expansions 

u{x,t) = f(a:,t)-^f"(x,t) + ^f^^\x,t)-... (5.1) 

vix,t) = _y/'(x,t) + |^/(3)(a,,t)-... (5.2) 

The effect of tlie surface tension in the free surface condition amounts to an 
extra term in 14.2|l : 



gh+^{u,+v,) + — ^x{t)^-^ (5.3) 



where T is the surface tension and x(i) again the arbitrary function depending 
only on time. 

This arbitrary function is ehminated by Boussinesq with the aid of the as- 
sumption that / and its derivations vanish for x — *■ ±cxd. Korteweg and De 
Vries drop this crucial assumption and x(t) is eliminated by simply differentia- 
ting (|5.3() with respect to x. It is now already noted that periodic waves are not 
a priori excluded from further discussion, this in contrast to Boussinesq, who 
used a different approach in his discussion of periodic waves, see section 7.1. 
Differentiating 1)5. 3|) with respect to x, the free surface condition 1)5. 3|l becomes 

g^j^u—+v— + ^^^--^-^ = {) (5 4) 

dx " dx " dx dtdx p dx^ 

Besides this we need again the kinematic condition H4.3|l 

dh dh ,^ ^. 

Korteweg and De Vries put ys = H + h{x, t) and f{x, t) = qq + /3(x, t) with 
go an as yet undetermined constant velocity. Substitution of H5.1|l and H5.2|) into 
1)5. 4(1 and H5.5|l gives in a first order approximation for h small and for a wave, 
progressing in the positive a;-direction, the expression 



h = h{x-{qo + VgH)t). 

Adding to the flow a velocity qo = —^/gH, we obtain the Lagrange steady 
wave solution with 

dt ' dt ' 

and 

9/3 qo dh g dh .9 /i, , ^ 

ox H ox qo ox qo 

where a is an undetermined constant. 

The next approximation is obtained by 

f{x, t) = qo - —{h{x, t)+a + ^{x, t)) 
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with 7 small in comparison with h and a. Substitution into (|5.4I) and (|5.5|l yields 
two equations for h and 7 and elimination of "/{x,t) gives finally the Korteweg- 
de Vries equation as it appeared for the first time in the thesis of De Vries |25) : 

dh 3 g d fl^^ , ^ I _d'^h 

with (T = - m, 

3 pg 

The addition of the velocity qo — a — —^fgH + j Ha to the flow may 
also be obtained by a transformation of the fixed (x, y) coordinate system to 
the moving frame 



^ = x-WgH- Jj^a% T = t. (5.7) 
Hence, in this moving frame and forgetting about the added velocity, we get 
dh i IT d [1,2 2 , 1 d^h\ „ 



dr 2\l H \2 3 3 ^ 

This equation with T = is equivalent with Boussinesq's result (|4.11l) and may 
be obtained by substititution of H5.7|l into (|4.11|1 . 



6 The Solitary Great Wave 

Because of the equivalence of the differential equations (|4.11|l and (|5.8() for the 
surface of a wave with amplitude, vanishing at infinity, it is evident that the 
theory of Boussinesq leads to the same results as that of Korteweg and De Vries, 
if capillarity is neglected. A necessary and sufficient condition for the existence 
of a steady wave is a constant uniform wave velocity in all points of the wave 
surface. 



6.1 The Solitary Steady Wave in the Theory of Boussinesq 

It follows from H4.8|l that the wave is stationary if w = \/gH + g/ Hhi, where 
hi is some as yet unknown constant, independent of x and t. Therefore, 



Integration with h ^ and ^ for x ±00 gives 



h{x, t) = /iisech^ 



3/11 



(x — ut) 



(6.1) 



(6.2) 



It follows hi > h{x, t) and hi is the crest of the wave. The wave velocity 

(6.3) 
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contains a correction of Lagrange's result with w = ^/ gH and was already 
experimentally verified in 1844 by Scott Russell The expressions H6.2|) and 
Hfci.3|l were also obtained by Rayleigh, see and H3.7|l . It appears that the 

wave velocity is the larger the higher the crest of the wave and this means that 
in the case of several separate solitary waves of the form ()6.2|l the higher waves 
will overtake the lower ones whenever the higher waves were initially behind 
the lower waves. It is known that this occurs without change of form, however, 
there is only a phase shift. 

The solitary waves behave like a row of rolling marbles, where the faster 
marbles carry over their impuls to the slower marbles. They were coined by 
Zabusky and Kruskal |3] as "solitons" to indicate their particle like properties. 
For an explicit calculation of this behaviour the interested reader may consult 
ref. 221' P^^'t II7 3.5. 

6.2 The Solitary Steady Wave in the Theory of Korteweg 
and De Vries 

Korteweg and De Vries do not have at their disposal an explicit expression for 
the wave velocity. However, for a steady wave in the moving (^, r) frame (|5.7|l 
one has §7 = and so by H5.8|l 

^U'^ +r^+rde^j^o (6.4) 

where a is the still unknown correction of the wave velocity. Integration under 
the assumption /i, ^ ^ for ^ ^ ±00 results into 



d/i _ -h^{h + 2a) 
di" V a ■ 

There are two distinct cases, ct > and cr < 0; we restrict our calculation to 
the case ct > since the other case can be treated similarly. If a > then the 
constant 2a is negative and taking 2a = — /12 one gets 



/^(O^tehMW-^e , (6.5) 




with /i2 > the crest of the wave. When the surface tension T is neglected, the 
parameter a reduces to a = ■^H^ and H6.5() is in agreement with Boussinesq's 
result (|6.2|l . Also the wave velocity is conform H6.3|l because we have by (|5.7|l 

uj = ^/gH - ,J^a = ^/gH +^,J^h2 (6.6) 
Korteweg and De Vries consider also solitary steady waves with a negative 



'3T 

C ^ U, l.K^. ±± ^ \! 

mately i cm for water. 



amplitude, possible for cr < 0, i.e. H < J ^ ; this quantity equals approxi- 
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7 Periodic Stationary Waves 
7.1 Boussinesq's Theory 

Boussinesq investigates in his Memoire "Essai sur la theorie des eaux cour- 
rantes" la Forme la plus generale des intumescences propagees le long d'un canal 
horizontal et rectangulaire, qui avancent sans se deformer" (^|, pp 390-396). 
By now, he does not have a formula for the wave velocity as in H4.8|l , because the 
condition h{x, t) — > for x — > ±oo does no longer hold in the case of periodic 
waves and so an easy evaluation of the wave form by setting lu constant is no 
longer possible. He introduces the mean horizontal velocity 



1 



H+h{x,t) 



which satisfies the relation 

see ^nii 276, pp 390-391, or the paper by St. Venant ,19, where a rather short 
derivation is presented. 

For a steady wave we have |^ = ^ = and ^ = -w^, so 

the latter equation becomes 

-ujU + g{h~H) + — + -^g^ = constant := ^c', (7.2) 

with uj the unknown constant wave velocity and ijjpc' the constant of integra- 
tion. 

The mean velocity U is eliminated with the aid of the conservation of flux 
in a reference system bound to the wave, i.e. 

H+h{x,t} 

udy — ujh{x, t) 

Jo 
or 

U = ^^^(l-h.), (7.3) 
H + h H \ h) ^ ' 

Substitution into H7.2|) and omitting terms of order and higher we get 
and by multiplication with ^ and integration with respect to x 



(^^\" = -^Ah' -Hh^(l- - dh - c"}, (7.4) 



\dx J \ LiP- 

with c" the integration constant. 
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Figure 3: Cnoidal Wave 

Suppose now that the minimum value ho of h, which is necessarily negative, 
is assumed at a: = and the wave surface represented by h{x,t) = h{x,t) — ho, 
then h > and ft, = at x = 0. This shift results into 



(7.5) 



where oq, ai and 02 are certain constants, which due to the condition h — and 
^ = at x = satisfy aq = and ai < 0. Therefore, H7.5() may be written as 

2 

^^{Ch + k)hil-h)}, (7.6) 

with k and / positive constants. 

Integration leads to a Jacobian elliptic function, but Boussinesq recommends 
to use Newton's binomium for a series expansion of the variable x into powers 
ofh. 

It follows from (|7.6|) that the solitary steady wave is only a particular case, 
resulting from (|7.6|l for k ^ 0. Finally, one may substitute h — h + ho in 1)7.4(1 . 
compare the coefficients of in H7.4|l and H7.6() and one obtains after a simple 
evaluation the following expression for the wave velocity 

uj^ ^g{H+{l-k)}. (7.7) 

7.2 The Theory of Korteweg and De Vries 

Since the Korteweg-dc Vries equation 1)5. 8|) . as derived in section 5, may also be 
applied to waves, not necessarily vanishing for x —^ ±00, the equation for the 
amplitude h{£^) of a steady wave is given by 

d /l.o 2 . 1 (Ph 
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Integrating this expression two times one obtains 



and 



1,, 2 , 1 d^h 

''^2^ +r''+rde=' 



C2 + 6cih + h'^ + 2ah^ + a i^— j =0, (7.9) 

with ci and C2 the constants of integration. 

The wave surface may be defined as y = Hq + h{^) with Hq the minimum 
value of y, h{^) > and h{0) = 0. 

It follows that ^ = and ^ > for /i = and so C2 = and ci < under 

the assumption a > 0. Consequently, the equation /i^ + 2afi + 6ci = has a 
positive root I and a negative root — fc and l|7.9|) reads 



^^±^lih + k)hil-h), (7.10) 

which for r = is the same as (17. 6|) . 

With the aid of the substitution h ~ I cos^ x Korteweg and De Vries obtain 
the well-known periodical "cnoidal" wave 



/i(0 = ;cnMJ^H, (7.11) 




where cn denotes one of the Jacobian elliptic functions with modulus M = , 
period 



AK^4 [ {l-t^y^il-- M'h^y^dt 
Jq 




and wave length AK^j^. 

This wave length becomes infinitely large for k —^ and the result is the 
solitary steady wave (|6.5(l . However, one gets for large values of k, i.e. for small 
values of M, the sinusoidal wave 



h{^) = I cos^ X — I cos^ 



with decreasing wave length for increasing k. This agrees with a result of Stokes 
[21] and in this case h{£,) may be expanded in a Fourier series; this may be the 
reason why Stokes at first believed that the only permanent wave should be of 
sinusoidal type. 

The approach of Korteweg and De Vries as given here is in particular attrac- 
tive, because of the central role of their equation H5.8() to which they frequently 
revert in the development of their theory. 
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8 The Stability of the Stationary SoUtary Wave 



It follows from the {x,t) dependence of the wave velocity ui, H4.8|l . that wave 
propagation involves in general a change of form, but by definition this does 
not occur in the case of a steady wave and so the question arises why the 
steady solitary wave is stable and an exception to the rule. For the possible 
existence of the steady wave a further investigation is required, in particular 
with regard to the "parameters" determining the stable behaviour. This has 
been carried out by Boussinesq and Korteweg-de Vries in quite different ways. 
The presence of the non-linear term and the dispersion term ^H^^^ is 
already an indication for a possible balance, furthering the stability of the wave. 

8.1 Stability in the Theory of Korteweg-de Vries 

The authors consider a wave form close to that of the steady solitary wave 

h{^) ^hsech^ip^) (8.1) 

where h and p are as yet arbitrary constants with p near (see (|6.5|) . 

The deformation of this wave is determined by the equation H5.8(l and sub- 
stitution of H8.1|l gives an equation for the evolution of the surface of the wave, 
given by y = /i(C,t): 

^ = 3^hp{Aap' - h) |-sech^(K) + l^^^} ^^MpO- 

(8.2) 

Taking a — iap^ ~ this equation becomes 

— = 3\/ -j-hp{Aap'^ - h)sech^{pO tanh^(K)- (8.3) 

UT \ H 

The choice p — and thus a — —^h results into ^ — and we get the 
steady wave H6.5|l . 

A numerical analysis of (|8.3|l shows that the wave in its course becomes 
steeper in front and less steep behind when p < and conversely when 

h 



This result is in contradiction with the assertion of among others Airy, that 
a progressive wave always gets steeper in front and less steep behind. This 
opinion is conceivable if the dispersion is neglected. 

8.2 Stabihty in the Theory of Boussinesq 

Boussinesq considers waves, not necessarily steady, with the same energy 



2 



-oo JO 



oo 



. ™ pH+h . ™ 

pgE = ^pg I h^dx + - dx {u^ + v'^)dy ^ pg h^dx (8.4) 
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see (El, pp 85-86). 

Furthermore, he introduces the functional 




which he calls the "moment de stabilite" and he shows that M is a conserved 
quantity, i.e. independent on t, ([^, pp 87-88, 97-99). 
After the transformation 

poo 

h^dx 



the expression for M becomes 

Boussinesq uses, without reference to Euler-Lagrange, the well-known method 
to obtain a condition for t) in order that M attains an extremal value and 
the result is 

2ij3. d f, dh^ 



From de = —h^dx — hda follows equation l|4.10(l with ^ — and therefore 
only the stationary solitary wave with given energy E yields an extremum for 
M. Variation of h with Ah gives M > for all h{e,t) and so the extremum of 
AI is a minimum. The stability of the wave is evident, because also AAI does 
not depend on t. 

As in discrete mechanical systems conserved quantities are of fundamental 
importance, also in continuous dynamical systems. Besides the integral invari- 
ants Q = hdx and E = h? dx, corresponding with the conservation 
of mass and energy, Boussinesq discovered a third invariant, the "moment de 
stabilite" . 

He also showed that the velocity of the centre of gravity of a wave does not 
depend on time and this implies a fourth invariant of the motion ([H], pp 83-84; 
[ISl, p 135). 

Conserved functionals may be considered as Hamiltonians in continuous dy- 
namical systems and they play there a role analogous to the Hamilton functions 
in discrete systems. These continuous dynamical systems have been investigated 
only rather recently; the first fundamental results have been established in the 
seventies by Lax 27^, Gardner 28 , Zacharov ,29, and Broer 'JD . Nowadays, 
there exists an extensive literature on this subject; a valuable introduction with 
many references is the textbook by P.J. Olver |31| . 

The Korteweg-de Vries equation is the prototype of an integrable system 
with an infinite number of conserved functionals, mutually in involution with 
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respect to a suitably defined Poisson bracket. In particular the Korteweg-de 
Vries equation may be represented as a Hamiltonian system in the form 



with 




dh'H is the variational derivative of the Hamilton functional Ti. and e a scale 
parameter, ref. [13 part I, ch2; part II, ch 5]. 

The first term in Ti. is the Hamilton functional for waves in the Lagrange 
approximation and the second term is the Boussinesq correction, given by the 
"moment de stabilite" M. 

Hamilton's theory for finite discrete systems dates from about 1835 and 
it was a century after Boussinesq that this theory has been generalized for 
continuous systems. Boussinesq has set, by using functionals, a first step into 
the direction of this generalization. 

9 Concluding Remarks 

We have discussed in the preceding sections the more important aspects of the 
work of Boussinesq and Korteweg-de Vries, who have besides these also studied 
other specific topics such as the velocity field, the path of the fluid particles and 
the motion of the centre of gravity of a solitary wave. 

Boussinesq finishes his article in the Journal de Mathematiques Pures et 
Appliquees with a qualitative examination of the change of form of long non- 
stationary waves and an attempt to prove that a positive solitary wave can be 
splitted into several other solitary waves. Korteweg and De Vries wanted to 
show that their approximation of the surface of a steady wave may be improved 
indefinitely, resulting in a convergent series. They claim in the introduction of 
their paper that "in a frictionless liquid there may exist absolutely stationary 
waves and that the form of their surface and the motion of the liquid below it may 
be expressed by means of rapidly convergent series". The calculations, however 
elementary, are so complicated and tedious that one may expect that these have 
not received much attention. Even the second approximation (p. 443), following 
on the formulae (|6.5II and (|7.11() of the present paper requires already so much 
effort that it is reasonable to be content with the first approximation as given 
in (jlSl) and f7TT|l . 

It is somewhat surprising that Korteweg and De Vries refer in their paper 
only to Boussinesq's short communication in the Comptes Rendus of 1871 [7| 
and not to the extensive article in the J. Math. Pures et Appl. |U] and the "Essai 
sur la theorie des eaux courantes" ^U] ™ 1872, respectively 1877. However, we 
should realize that the international exchange of scientific achievements in those 
days was not at the level as it is today. 
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As to the credit of the "a priori demonstration a posteriori" of the stable 
sohtary wave, this credit belongs, of course, to M. Boussinesq. On the other 
hand, Korteweg and De Vries merit to be acknowledged for removing doubts on 
the existence of the "Great Wave" and for their contribution to the theory of 
long waves in shallow water. 
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